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I. 


INTRODUCTION 


When  a distant  target  is  illuminated  by  a coherent  light 
source,  the  far-field  backscatter  intensity  is  speckled.  In  general, 
there  is  a constant  component  and  a spatially  varying  component. 

Hence,  the  intensity  measured  by  a backscatter  detector  (at  the 
illumination  site)  fluctuates  as  the  target  body  turns.  If  the  target 
is  a flat  surface,  a discrete  Fourier  transform  (DFT)  of  the  back- 
scatter intensity  indicates  the  decorrelation  length  for  the  random 
phase  variations  of  the  reflected  light  at  the  target  surface  (due 
to  surface  roughness,  for  example  [1,  2]).  If,  however,  the  body  is 
lurved  in  such  a way  that  the  surface  always  moves  through  the  same 
surface-of-curvature  (e.g.,  a cone,  cylinder,  or  sphere  rotating 
about  the  axi s-of-symmetry) , then  the  DFT  will  only  yield  information 
concerning  the  target  width,  provided  its  surface  is  rough  [3]. 

It  is  for  this  class  of  targets  that  the  present  formulation  is 
developed,  and  it  is  pertinent  to  monostatic  laser  radar  applications. 
The  ratio  of  rms  Intensity  fluctuation  to  average  intensity  is  formed 
from  the  measured  data.  This  ratio  is  dependent  upon  the  rms  random 
phase  fluctuation  of  reflected  light  at  the  target  surface  and  the 
number  of  decorrelation  cells  on  that  surface.  For  the  case  where 
the  phase  fluctuation  is  due  to  surface  roughness,  the  equivalent 
parameters  are:  (1)  rms  height  fluctuation  and  (2)  the  product  of 

average  roughness  slope  and  the  square  root  of  the  exposed  target 
surface.  Both  of  these  parameters  can  be  uniquely  determined  if 
measurements  are  made  as  a function  illumination  wavelength. 

Figure  1 illustrates  the  measurement  configuration.  The  angle  a 

is  the  orientation  of  the  turning  surface  at  any  instant.  In  the 

model  presented  in  this  treatment,  the  target's  surface  field  is 

assumed  to  vary  from  the  average  in  a half-cycle  sinusoid  within  each 

decorrelation  area  cell.  The  amplitude  and  period  for  each  cell  is 

independent  of  that  in  any  other.  The  effective  phase  for  each  cell 

Is  random  with  equal  probability  between  the  bounds  - ({i  and  <p  , 

m m 

The  phase  fluctuation  statistics  are  assumed  to  be  the  same  for  all 

2 

directions  over  the  hody.  The  decorrelation  area  As  = il  defined 

c c 

here  relates  only  to  fluctuations  of  the  reflected  surface  field  from 
the  average.  It  is  further  assumed  that  the  number  of  decorrelation 
cells  on  the  body  Is  large. 

Polarization  is  not  explicitly  treated  in  this  report.  If  the 
polarization  effect  relates  to  the  surface  curvature  (from  the  line  of 
observation),  forming  the  ratio  of  rms  intensity  fluctuation  to 
average  Intensity  cancels  out  its  effect. 


The  Doppler  effect  due  to  target  rotation  is  not  included  in  this 
treatment.  Although  the  far-field  intensity  due  to  a given  point  on 


the  surface  Is  critically  dependent  upon  the  Doppler  shift,  there  is 
a compensating  effect  in  the  integration  of  the  effect  due  to  the 
whole  surface  which  greatly  reduces  its  effect  on  tlie  overall  statistics 
of  the  far-field.  Thus,  the  treatment  presented  here  is  valuable 
despite  the  exclusion  of  the  Doppler  effect.  This,  however,  does  not 
imply  that  the  Doppler  effect  is  not  important  to  the  spread  in  the 
frequency  spectrum  of  the  reflected  light. 

The  present  work  is  intended  to  be  a guide  to  the  understanding 
and  development  of  laser  radar  applications. 


II.  THEORY 

A.  General  Formulation 

A target  body  which  rotates  so  that  its  surface  always 
passes  through  the  same  surface-of-curvature  is  considered.  Examples 
may  be  a snhere,  cone,  or  cylinder  rotating  about  their  axes  of 
symmetry.  The  axis  is  assumed  to  be  approximately  perpendicular  to 
the  line  of  observation.  The  distant  target  is  illuminated  with 
coherent  light,  and  the  far-field  backscatter  is  detected  at  the 
illumination  site.  The  electric  field  amplitude  at  the  detector  is 

A = A + AA  , (1) 

a 

where  A^  is  the  average  amplitude  and  AA  is  the  fluctuation  from  the 

average  which  occurs  as  the  target  rotates.  Thus,  a speckle  field  is 
constantly  sweeping  across  the  detector.  According  to  the  Rayleigh- 
Sommerfeld  diffraction  formulation  [4],  A is  related  to  the  surface 
field  a of  the  target  by 

A = Ja  exp(jk(r  - R)]  cos(n,£)  ds  , (2) 

S 

where  R is  the  distance  from  the  detector  to  the  nearest  point  of  the 
target  surface,  r is  the  distance  to  the  surface  at  any  given  point 
(average  surface  position  at  this  site  if  the  surface  is  rough) , and 
cos(n,r)  is  the  angle  between  the  surface  normal  and  the  line  of 
observation  (Figure  1).  The  integration  is  over  the  exposed  surface. 
It  is  assumed  that  the  integrand  of  Equation  (2)  can  be  expressed  as 

(r)  F (n , r ) , 


6 


where  is  cos(n,£)  times  a phase  shift  factor  related  to  surface 

curvature  (a  phase  component  which  is  due  to  the  fact  that  different 
parts  of  the  target  surface  are  not  equidistant  from  the  plane  of  the 
incident  illumination  field).  The  is  the  field  amplitude  which 

would  exist  at  the  surface  if  the  curvature  could  be  neglected  and  thus 
only  contains  a phase  factor  determined  by  the  detailed  surface 
properties  such  as  roughness,  dipole  irregularities,  etc.  The  a^ 

may  depend  on  the  angle  n^,r^  in  some  cases.  For  example,  if  fluctua- 
tions in  are  due  to  roughness,  then  an  exponent  2k  h(r)  cos(n,£) 

is  Involved  (h  = height  from  average  surface).  However,  there  are 
compensating  aspects  that  allow  one  to  assume  an  effective  phase 
change  due  to  roughness  of  2k  h^^  for  the  i th  decorrelation  cell 

(Appendix  A).  The  assumption  is  made  that  the  phase  of  the  cells  does 

not  depend  on  the  surface  orientation  and  fluctuates  randomly  with 

equal  probability  between  <fi  and  - . Thus, 

m m 


= a^  p expjif'  , (3) 

where  a^  is  a real  constant  representing  the  plane  wave  amplitude  of 

the  incident  field,  p is  the  reflectivity,  and  ifi  Integrated  over  a 

decorrelation  cell  varies  randomly  between  ({>  and  - <}>  , depending 

m m 

only  on  the  region  of  the  surface  being  considered  (not  its  angular 
orientation).  Thus,  is  the  surface  field  on  the  surface  if  it  were 

flat  and  normally  Illuminated. 

The  terms  of  Equation  (1)  are  given  by 

and 

AA  = expjkR  j Aa^^  F ds  , (5) 


where 


la 


<ai> 


(6) 
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and 


Aai  = a^ 

- a 

la 

(7) 

The  quantity 

AQj^  fluctuates  as  a random  function  of  a. 

It  is  seen  that 

/-la 

ds  = J F ds  , 

(8) 

since  a,  is  a constant  over  the  surface, 
la 

The  decorrelation  cell  is  defined  as  the  constructed  surface  area 

for  which  a,  locally  departs  from  a.  , either  positively  or  negatively 
X 

(Figure  2).  Thus,  the  decorrelation  cell  area  is  not  necessarily 
proportional  to  the  square  of  the  correlation  length  based  on  an  auto- 
correlation of  ij-  Such  irregularity  in  the  behavior  of  41(01)  exists 
that  a correlation  of  Aa^  with  itself  will  fall  off  to  an  average  of 

zero  for  displacement  of  Even  if  Aa^^  is  varying  sinusoidally 

[with  higher  frequency  than  4>(H)1  because  X ''  k 4i  (4i  = 4>  ).  the 

m m max 

derivative  d4>/di  is  assumed  to  vary  sufficiently  that  the  period  of  the 
sinusoid  fluctuates  from  one-half  cycle  of  Aaj^  to  the  next  so  that  the 

correlation  function  reaches  its  first  zero  at  about  !!,=  )!,  and  then 

c 

oscillates  randomly  with  diminishing  amplitude.  Thus,  each  "peak" 
or  "valley"  of  Aa^^  is  held  to  be  independent  of  any  other. 

It  is  important  to  realize  that  the  decorrelation  area  defined 
here  refers  to  random  fluctuations  in  Aa^^,  however  small  compared  to 

and  not  phase  or  a^^  Itself.  A reordering  of  the  Integration  in 

Equation  (5)  is  performed  such  that  Aa^^  runs  smoothly  from  the  highest 

to  the  lowest  value,  and  F fluctuates  randomly  about  its  average.  The 
number  of  cells  N is  taken  to  be  sufficiently  large  that  the  fluctuating 
F integrated  over  an  incremental  range  of  (small  compared  to  the 

maximum  value  of  Aa^)  is  the  same  as  the  overall  average.  This  allows 
one  to  write* 


To  be  absolutely  rigorous,  the  Integral  over  Aa^  is  not  independent 

of  F.  However,  the  "average  rms  fluctuation"  of  E ds/<F>  as  a 

varies  is  the  same  as  that  for  the  integral  over  Aa^  alone.  Since 

tills  treatment  leads  toward  measurements  of  averages.  Equation  (9) 
may  be  assumed  without  error. 
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(9) 


F ds  ~ <F>J 


ds 


where 


= U' 


<F-  = -|  / F ds 


(10) 


S is  the  total  exposed  area  of  the  target.  A clarification  of  Equation 
(9)  is  shown  in  Figure  3. 

Note  that 


ds  - ^ Aa^.  As.  , 


(11) 


i=l 


where  is  the  average  value  of  Aa^  in  the  i th  member  of  the  N cells 

and  As.  is  the  correlation  area  of  that  cell.  The  area  As.  fluctuates 
1 1 

independently  of  about  an  average  value  As^.  Hence, 


/Aai 


ds  = As  5Aa-  = As 
cl  c 


S ^^li 


i=l 


or,  using  Equations  (9)  and  (10), 

As 

ds  6A3j^ 


Jf  Aa^  ds  - J F 


(12) 


(13) 


where 


6Aa 


1 1 ^^i  ’ 


1 = 1 


and 


(14) 


As  /S  = N 
c 


(15) 


Thus,  using  Equations  (4),  (5),  (8),  (9),  (10),  and  (13), 


A . / F ds 

a lAR  la 


/• 


(16) 


11 


:egratio! 


and 


AA  = i 6Aa,  f F ds  . (17) 

3AR  N 1 J 

Therefore,  I = g A*A  is  given  by 
or 

■ « (n/  I /■’  [“L  * f “la  ‘“1  * 7 

where  the  fapt  that  6A3j^  + AAa^^*  = 2 AReAa^ . I is  the  far-field 

intensity  which  fluctuates  as  the  target  rotates.  Ttius,  one  may  write 

1(a)  = ^1>  + 6Ij^(a)  + 6X2(01)  , (20) 


where 


= I/f  dsl^  |a2^  <|6AaJ^>j 

2 

'h  = 8 (j^)  1 /f  dsl^  I a^^  6ReAa^  , 


(21) 

(22) 


and 

“2  - « {rj  I ^ “"I'  ~2 

Note  that,  since  it  fluctuates  about  zero,  <6ReAaj^>  = 0.  The 
averaging  relates  to  the  variable  a. 

2 

Now,  turn  to  the  derivation  of  6ReAa^  and  |6Aaj^(  in  terms  of 

surface  properties.  First,  consider  that  the  fluctuation  Aa^^  in  the 
n th  cell  is 


Thus , 


6Aa,  = > Aa,  = a p / (expj<J>  - 'uxpj  ij'>) 

i ^ in  o ^ n 

n=l  n=l 

Now,  choose  the  set  such  that  fiAa^  ' ~ h,  and  ri-order  tiu' 

summation  to  produce  the  set  {41.}  which  runs  from  liiglii-st  to  lowest, 

but  includes  the  same  values  as  {<J>  }: 

n 


6Aajy  = 0 = a 


p > (expj({)  - <expj4)>) 

O i. 


Now,  let  the  surface  orientation  a change  so  tliat  4 • 4).  + by  ( .A.ij 

is  no  longer  zero).  Thus, 


6Aa,  = a 
1 o 


P ^ ^expj(4>^  + - <expj(J)>J 


Reordering  the  terms  such  that  the  summation  runs  from  higliesl  14^  + 

to  lowest  (but  keeping  the  same  values)  and  calling  tliis  set  H||',  iSAa^ 
is  expressed  as 


>1  = Ho  P ^ (expj())^  - <expj< 


>) 


Due  to  the  reordering,  the  absolute  difference  between  41^^  and  4^  (i  * h) 
is  I'54'j^I  which  is  much  smaller,  on  the  average,  than  |'54'.|-  Tiieretore, 

N 

6A3j^  = SAa^^  - SAaj^g  ~ ^ [^expj(4i^  + ^<l>|j)  - t‘xpj4  J 


i=h=l 


becomes 


p ^ ^expj4’j  (i  + ” ^expJ4>>J 


i=h=l 


• l “ p ^ (expj4'^) 


i=h=l 


Mso  has  been  made  of  the  fact  that  the  sum  over  expj<{i.  is  N <expJ4>>. 
Thus , 


N 


Im  6Aa 


, = a p > cos<t>. 
1 o 1 


641, 


i=h=l 


Im  6Aa,  = a p 
1 o 


^ cos4>.  - ()).)^ 


'^h  ^i>os 


pos 


S cos^.  (4^  - 


neg  4 


N/2 


Im  6Aa,  = a p <cosd!>  > (4),  + 4>  ) , 

1 o ^ h pos  n neg 

h=l 

where  the  subscripts  pos  and  neg  denote  the  value  for  <|ij^  for  the  h th 
cell  (positive  i()^  region)  and  the  h + N/2  th  cell  (negative  41^  region), 
respectively.  Defining 
N/2 


(4>,  ■*'4),  ) 

h pos  h neg 


(27) 


h=l 

enables  one  to  write 


Im  6Aa,  = a p <cos4>> 

1 o I 


(28) 


Also,  since  E641,  = 24>,  > 

h h 


V = S 


j-1 


where  {4j}  Is  a random  reordering  of  { 4>,  ) • According  to  the  formula 
J h 

for  the  fluctuation  of  a sum  of  N fluctuating  terms,  the  rms  value  of 
6j4>  over  all  orientations  a is 


(6,4>)  = N^^^  (4>.) 

I rms  j rms 
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Due  to  equal  probability  of  all  ({).  ^ , the  size-ordered  array  is 

linear,  and  it  is  simple  to  show  that  (<P . ) =4'  />^  , where  I is 

' rms  m 


the  maximum  value  of  <})j  in  the  set  of  N(u)  values  whieli  lias  been 
averaged  over  all  orientations  (all  values  of  a).  Thus, 


(6  A)  = <p  //T 

I rms  m 


Attention  is  now  turned  to  Re6Aaj^.  From  Equation  (26) 


Re  6Aa  = _ a p \ siniji 
■t  o i 


l=h=l 


Re  6Aa^^  = - a p <sln|(()l>  fl  ‘♦h  - I 64>j 


pos  neg  4)^ 


Re  6Aa^  = 


a^  p <sln 141^ 


^ ^"^h  pos 


pos  4 
i=h 


neg 


neg  4 
i=h 


Re  6Aa^  s - a^  p <sln|4|>  ^ ('P 


k - 4,  ) - N < 4. 

pos  h neg  ' i 


Thus,  defining  <5  4 as  the  term  in  brackets, 
K 


Re  6Aa  » - a p <sin  <j)|>  6 $ 

i O iv 
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Now,  define  a set  iji!  such  that  (fi'  = <|i,  for  i = h.  Equation  (32) 
1 in 

becomes 

N 

Re  6&aj^  = -a^  p <sin|4>j>  ^ (U[|  - U^|) 

i=l 

Since 

N 

2 = N <U|> 

i=l 

then 

N 

2 u;i  - N<ui> 

i»l 


2 

By  forming  <6  > and  noting  that,  for  the  rectangular  probability 
R 

function,  <I<1>|>  = *^/2  <i>  , it  is  straightforward  to  show 

' ' rms 


= ♦ //l2  = i (6^0 

m z I 


rms 


(34) 


Wliereas  6^;})  is  a measure  of  whether  a particular  set  } inclines 
1 n 

toward  positive  or  negative  values,  6 A is  a measure  of  whether  the 

R 

total  range  of  is  greater  or  less  than  that  of  A.  study  of 

Equations  (27)  and  (32)  makes  this  clear.  Thus,  6 4>  and  6 ())  both 

1 K 

fluctuate  about  zero  and  have  similar  rms  value,  but  they  are 
Independent  (random  with  respect  to  one  another) . 

I 1 2 

Therefore,  | 6Aa^ | is  given  by 

|6AaJ^  = a^  p^  T cos<> (6^.^)^  + <sin|4-|>^  1 , (35) 


and,  hence,  defining 

® « (n)  ‘*"1^ 


(36) 
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and  noting  that  a,  = > = <cos())>  a p.  Equations  (22),  (21),  (11), 

I *1  I o 

(35),  and  (36)  give 


(6I)j  = 


2 2 1 

-B  2 a^  p — <cos4)><sin  1 ())  I > 


(17) 


and 


(61)2  = B a 


2 2 / 1 f 
o p u; 


, <cosd)> 


+ <sin|(|)l>^ 


j^(6j4>)^  - ' (6j<t>)^>j 

,*)^>]  j • 


Because  6I  = (61)^  + (61)2, 


(61) 


rms 


+ 2 (6I)^(6i)2  + 


(‘0^] 


1/2 


(18) 


(39) 


The  product  (61)  (6I)_  vanishes  upon  taking  the  average  since  6 <p 

it  R 

fluctuates  about  zero  and  (6I)„  is  a linear  combination  of  constants, 
2 2 

and  (5  <())  , and  (6  (ji)  which  is  independent  of  6„(j).  Thus, 

K 1 R 


1/2 


(40) 


2 2 2 

To  use  Equation  (40),  the  expression  <l(6i}i)  - <(6ifi)  >]  > , where  6 
may  be  subscripted  with  1 or  R,  must  be  reduced  to  one  with  known 
parameters.  Note  that 


<[(6<fr)^  - <(6(<.)^>]^>  = <(S<t)^  - 2 (6<p)^  <(6<t>)^>  -h  <(S(P)^>  > 


4 2 

= <(641)  > - <(64>)  > 


Appendix  B shows  that  for  6^  • 6||r((i,  this  further  reduces  to 


(N‘ 


N) 


2 


(41) 
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•9> 


T 


Since  <<()?>  = = (<l>  />^)^  = 4i^/9  for  6 <(i,  and 

1 1 rms  m tn  1 

2 2 

<(6j^)^>  - <(6j.(^)^>  = 16<(6^(^)^>  - <(6j^«)^>  = (N^  - N)  ^^/9  . 

(42) 

Hence,  using  Equations  (34),  (37),  (38),  and  (42), 

2 

<(6I)^>  = 4 a^  <cos(ji>^<sin I (fi  I N (j>^/12  (43) 

1 o \N/  m 


and 

<(6I)2>  = a^  ^(<cos(fi>^  + <sin|(J>|>^)  (N^  - N)  4'^/9j 

(44) 

The  cross-term  2/16  <cos(fi>^<sin  | 4i  | >^  [(6j({i)^  - < (6ji}i)^>]  [ (6j^4))^  - <(6j^(J))^>] 

has  been  omitted  because  each  factor  fluctuates  about  zero  and  is 
independent;  hence,  it  vanishes  upon  averaging. 

2 2 2 2 

Noting  that  = a^  p <cosij)>  and  using  Equations  (21),  (34), 

(35),  and  (36),  <I>  is  given  by 


<I>  = B 


2 2 
p <cos<t>>  + 


(If 


a^  p^  (<cos<}i>^  + ^ <sin|<Ji|>^)  N i)>^/3 


or 

<I>  » B a^  p“  ^<cos$>^  + ^ (<cos<{)>^  + <sin|$|>^)  . (45) 

Using  Equations  (40),  (43),  and  (44),  (61)  is  given  by 

rms 

^^^^rms  " ® ‘^m^^  |<cos^>^<sin U | >^ 

2 

+ (<cos4i>^  + <sin|4>|>^)((|>^/3)  (N  - 


or 
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(61)  = B 

rms  o 


|<cos(^i>^<sln  I ^ I 


« o a p TiZ  pcus(p-' 

° (3N)^''^  I 

+ (<cos.^>'‘  + ^ <sinU|>S  (0^/3)  (i  - 


1/2 


. (46) 


Now, 


<COSl}l> 


cos4> 


d$ 


sinit 

m 


and 


<s  in  I (j)  I > 


d<p 


1 - cosd) 

m 


The  spread  of  i is  + i}i  //fT  , thus,  for  N > 100,  this  spread  is 
, . . , , ra  m 

negligible. 


Thus,  from  Equations  (A5)  and  (A6) , 


(A7) 


(A8) 


<co8*>^<ElnU  1 + (<cos^>‘*  + <sin|^|>^)  (n  " ~t] 

r 

<co8*>^  + ^ (<cos^>^  + 

(49) 


Thus,  the  ratio  of  far-field  intensity  fluctuation  to  the  average 
is  uniquely  determined  by  the  rms  phase  fluctuation  = the 

number  of  decorrelation  cells  N = S/As^  . Appendix  C presents  a 

derivation  of  N in  terms  of  surface  properties. 

Case  Where  is  Due  to  Surface  Roughness;  A scheme  will  now  be 

developed  to  deduce  the  rms  height  variation  and  average  absolute  slope 

of  a surface  (where  (p  = -2kh)by  measurement  of  at 

m m 

different  wavelengths.  First,  Imagine  a sweeping  through  a broad  wave- 
length band.  Whenever  <p  = ^ v h /X  * it,  then  <cos!|)>  ■>  0 and 
“ mm 

Equation  (49)  becomes 


1 

|<sinl<ti|>^^ 

- - — 

(61)  „l/2 

rms  . N 

3 

nvJ 

■ - ■ ->  

<I> 

1 . ( 1 2 

- <sin  ((«>(> 

3 
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which,  for  large  N,  is 


V 1 X > X = 4 h = 4 ^ h 

<1>  ’ 1 m ris 


Thus,  X,  = 4v^h  is  the  condition  on  X for  which  first 

1 rms  rms 

reaches  unity  as  X is  decreased  from  larger  values.  This  condition 

uniquely  determines  h 

^ ' rms 

Now,  imagine  X is  increased  well  beyond  2 X^^.  In  this  region, 

A « n/2  (or  h «X_/8);  hence,  1 is  given  by  Equation  (C-6)  of 
m mi  c 

Appendix  C.  Thus, 


N = S/As  = ^ S < ^ > /h^ 

c c dil  m 


Also,  due  to  smallness  of  <}>  and  1/N,  Equation  (49)  becomes 

m 


r<cos(|>>^<sin I aI  >^1 


<COS(J)> 


rms  _ m 


1 <sin|4>l> 

1/2  <cos())> 


or.  since  A is  small  and  (|>  = 4nh/X-, 

’ mm3 


4 Tt  h 


X3  (3N) 


1/2  \ 1 


4 ir  h 


Further,  using  h = /T  h and  Equation  (51), 
^ m rms 


3(4tt)  / rms 
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An  intermediate  value  of  A,  e.g.,  A^,  is  chosen  such  that 

(t  = tt/2.  At  this  point,  A„  = 8 h = S/T  li  and  ?;quation  (C-6) 
^m  ‘2m  rms 

again  yields  Equation  (51).  Also, 

<cos(}i>  = <sin|ijil>  = 2/tt 


Hence,  for  large  N,  Equation  (49)  becomes 


or 


(53) 


Any  two  of  Equations  (50),  (52),  or  (53)  can  be  used  to  solve  for 

h and  <|dh/dS.l>: 
rms  ' 


(54) 


and 


also. 


idhi. 

2 

//(6I)  \ 

7T 

2 . 

I rms  j 

64/T 

\ /at  A3 

(55) 


(56) 


where  = A^/2  is  the  largest  value  of  X for  which  ~ 

and  A^  is  any  value  of  A which  Is  much  greater  than  A2* 

A sketch  of  the  behavior  of  (4l)  /<!''  with  A is  shown  in 

rms 

Figure  4.  Calculations  used  to  aid  in  this  illustration  are  presented 
in  Appendix  D. 
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Figure  4,  Behavior  of  T •=  (61)  /<!>  versus  NOTE:  The  per- 

® rms 

forated  line  is  the  behavior  of  P when  is  replaced  by  a spread 
so  that  the  probability  density  function  for  <t>^  does  not  have  a 
sharp  cutoff  at  <t>  . In  the  figure, 

^2  = ^ <|dh/di|>) 

and 

r = h^„  1 + InCSyTh  /\  )r/(4S  <|dh/d^|>^)  . 

o rms  I ' rms  o I ' 

N is  assumed  to  be  very  large. 


m. 


DISCUSSION  AND  CONCLUSION 


The  far-field  backseat  ter  intensity  measurements  for  a 
particular  class  of  rotating  targets  illuminated  c<jl)erent]y  from  the 
detection  site  are  indicative  of  certain  surface  detail  features. 

This  class  includes  cones,  cylinders,  and  spheres.  If  ^ 

is  measured  as  a function  c>f  wavelength,  then  the  surface  rougliness 

) /2  I , 

4>  = k h can  be  obtained.  Also,  S “<  dh/di  > can  be  obtained, 

rms  rms  ' ' 

It  is  not  possible  to  determine  the  average  roughness  slope  unless 

the  exposed  target  surface  area  S is  known.  This  may  be  estimated 

from  radar  cross-section  measurements  or  Fourier  analysis  of  the  far- 

field  intensity  (which  only  yields  overall  size  information  for  this 

class  of  targets).  For  a cone  or  cylinder,  S = ti/2  S , wliere  S is 

3 3 

the  cross-section  area;  and  for  a sphere,  S = 2 S . Thu.s, 

S = 1.7  holds  for  a variety  of  sltapes. 


Provided  that  the 
terminates  abruptly  at 


distribution  function  for  (f . is  constant 

2 ^ 

• 1>  , then,  for  large  N (({i  /3N  » 1),  the 
m m 


but 


oscillations  of  f as  shown  in  Figure  4 for  A < A are  valid.  However, 

— o 

it  is  unlikely  that  the  distribution  function  abruptly  drops  to  zero 

for  I (|i  I > >1)  . In  practice,  the  value  (}i  should  be  replaced  by  a 
m m 

weighted  spread  of  ij)  about  <))  . For  a sufficiently  large  spread, 

m 

<cosii>  remains  approximately  sin'))  / until  <j)  >■  tr . Beyond  this 

mm  m 

point,  <cos4i>  is  very  nearly  zero.  However,  ■'sin|())|'>  maintains  a 
larger  value  than  <cos((i>  and  no  longer  has  periodic  zeros.  The 
result  is  that  the  expression  for  f = very  much  the 

same  as  derived  except  that  f ~ 1 for  all  A < 


If  the  reflectivity  p is  not  a constant  but  varies 
0 (not  random),  then  it  may  be  "lumped"  with  F under  tiie 
integral.  Since  the  surface  integral  over  F cancels  out 
the  ratio,  it  is  not  important  in  this  analysis. 


smoothly  with 
surface 
in  forming 
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Appendix  A.  JUSTIFICATION  THAT  0^  -2  k <h>  OVER  A 
DECORRELATION  AREA  CELL  WHEN  0 ^0  AND  PHASE 
IS  DETERMINED  BY  SURFACE  ROUGHNESS 


Consider  a constructed  surface  [perforated  line  in  Figure  A-l(a)] 

whicli  represents  the  average  surface  observed  at  an  angle  0.  Examine 

one  decorrelation  cell  for  X > 8 h (where  the  surface  departs  from 

m 

h = 0 and  then  returns).  Reflection  from  the  real  surface  at  distance 

r appears  to  come  from  the  constructed  surface  at  distance  r.  The 
r 

path  difference  between  the  wave  from  the  real  surface  and  the  wave 
if  it  struck  the  constructed  surface  is 


Ap  = r - r - H 
r 


Now, 


- h 


- h 


sin 


(e  - 4-J  sin(e  - 20.  - 


and 


r - r = H sin  i.  = 
r 1 


— ■ --  - ■ sin ('20  . - 

sin^0  - 20.  - ^ ^ 


Thus , 


Ap 


s in  ^0  - 20 
or,  using  trigonometric  identities. 


7Jsln(20i  - - 1] 

i 2/^ 


2 

cos  0 , 


Ap  = - 2 h 


cos(20^  - 0) 


This  does  not  depart  significantly  from  - 2 h except  for  0 2 ti/4, 
in  which  case  (for  small  roughness  slope  relative  to  0 so  that  0 = 0), 


'Ap>  = - 2 cos0^ 


(A-1) 
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tvnamfm 


The  Ray leigh-Sonunerf ield  Diffraction  formulation  yields 


A = J a [exp(jkr)/r]  cos9  ds  , 


where  the  surface  field  a is  given  by 


a = a^  p exp[jkr  + j<j)(r) ; 


(A-2) 


(A- 3) 


where  a^  is  a real  constant  representing  the  incident  field,  p is  the 

reflectance  (Includes  a factor  = - 1 if  the  surface  is  metallic 

and  inverts  the  phase  by  180  deg),  and  0 is  the  phase  of  a at  the  con- 
structed surface  due  to  roughness.  Hence,  dividing  the  surface  into  N 
decorrelation  cells  and  using  Equations  (A-1)  and  (A-3) , 

N 

A = ^ As^  [exp(jk<r>^)/r]a^  p exp[jk<r>^  - j2k<h>^  cosS^l  cosO^ 

i=l 

Examination  of  Figure  A-l(b)  will  show  that  the  intensity  of  the 
reflected  light  at  the  constructed  surface  peaks  above  the  projection 
of  the  cell's  midpoint  on  the  constructed  surface.  Thus,  is  a 

weighted  average  and  is  equal  to  r^^  + Ar^^,  where  r^  is  r at  the  projectec. 

midpoint  on  the  surface  and  Ar^  is  a small  value.  If  h is  negative, 

Ar^  is  positive;  and  if  h is  positive,  the  intensity  peaks  below  the 

projected  midpoint  and  Ar^  is  negative.  The  magnitude  of  Ar^  increases 

as  h increases  and/or  cos©  decreases.  Thus,  A becomes 


A 


^o  P expJkR 
j A R 


I 


As  . 
1 


exp[ j2(kr^  - R) ] 


X expj  (-2k<h>jj^cos9^  + 2k  Ar^)  cos9^ 


(A-4) 


The  effective  cell  phase  is  defined  as 


2 k cos6^  + Ar^(0j^)J 


(A-5) 
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IP 


Duo  to 

as  fast 

cient  ly 
thon 


tlie  compensating  nature  of  Ar.,  does  not  fall  off  with  cos0 
as  2k  •:h>  cos9.  The  assumption  is  made  that  for  cos0j^  suffi- 
large  that  the  contribution  to  the  integral  is  significant, 


(j) . ~ - 2k  <h> . 

i 1 


(A-6) 


Thus , 


A 


a p 


I 


As . 


where 


(A- 7) 


F(r^,0i)  = exp[j2(kr^  - R) ] cos0^ 


(A-8) 


is  defined  as  the  surface  curvature  factor. 

It  is  also  assumed  that  the  cell  slope  deviation  from  0 Is  suffi- 
ciently small  that  shading  due  to  height  fluctuations  is  not  signifi- 
cant for  COS0  sufficiently  large  that  contribution  to  the  surface 
integral  is  important. 
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below.  (An  observer  at  0 detect 


Appendix  B.  EVALUATION  OF  <(60)^  > - <(60)2  > 

Let  ()>^  be  the  phase  of  site  1.  Also  rake  6(p  as  representing  6^.1; 


60  = ^ 


i = l 


Thus , 


N 


N 


C?)^  = ^ 0^  + ^ 1 ^ 
i=l  i=l  i^L 


(B-1) 


(B-2) 


N 

(60)^  = N<0^>+^  \ <j>.  0.  , 

I ^ ^ 1 j 

i=l 

so  Chat 


(B-3) 


2 2 
<(60)  > = N <0^> 


Also , 


(60) 


Thus , 

C3i 


<(60)  > = 


N ^4*.^  + 0 + <f  ^ ^"1  ^ 


(B-4) 


J (l,\\ 

i=i  m \i=i  / \j^i  7 


(B-5) 


where  the  factors  of  the  last  term  average  separately  because  0^  and  0^ ^ 

(j^l)  are  independent.  Thus,  using  Equations  (B-1)  and  (B-4),  the 
last  equation  becomes 
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(B-6) 


(B-7) 


Appendix  C.  AN  EXPRESSION  FOR  N = S/As^ 


Tiu’  quantity  ^ is  the  average  of  the  maximum  value  of  <p . in  all 
m2  t 

decorrelation  area  cells  As.  piecewise  linear  so  that 

the  surface  field  Re  a^^  varies  sinusoidally  over  a cell  as  sliown  in 

Kiguri'  C- 1 . riTe  amplitude  and  periodicity  is  random  from  cell  to  cell. 

For  small  the  periodicity  of  Re  a,  is  the  same  as  that  of  ())  .. 

mi  '1  mi 

For  large  Re  a^^  may  oscillate  through  many  cycles  as  <p ^ changes 

and,  hence,  i . is  small. 

Cl 

Since  4>,  fluctuates  randomly  about  zero  from  - d>  to  (b  , then  i()  , 
i m m ml 

fluctuates  from  0 to  2 . 

m 

4 = (1  + 6!)  4,  , -1  < s!  < 1 , (C-i) 

ml  1 m 1 


where  5!  is  a bounded,  random  fluctuation  from  cell  to  cell,  and  the 
1 

average  cell  length  is 


1 


-1 


(C-2) 


Defining  t, 
0 to 

Re  Aa^^  Co 

small  <()  ). 
m 


as  the  surface  distance  for  which  f rises  linearly  from 

then  one  can  see  that  8,  . is  the  distance  from  "peak" 

pvi 

the  "valley"  of  Re  Aa^  in  the  region  of  the  i th  cell  (for 
Therefore, 


Small 

4>  condition: 

0 < , 

= (l  + 6')((>  <TI=> 

If 

m 

ml 

i m 

cl  pvi 

Large 

4>  . condition: 
ml 

• " Vi 

• (l  + 6p(()m<<»=> 

Restating  these  conditions  along  with  that  of  Equation  (C-1)  gives 


i) 

4>  < 

ir/2  condition: 

D.  = Z . ; 

m 

ci  pvi 

il) 

4>ni  > 

m 

n/2  , small 

condition:  -l<'5l<';r !■■>  , 

1 4 cl  pvi 

m 
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ill)  d)  > ir/2  , large  6l  condition:  - 1 < iS ! < 1 =>  It  . 

m 1 Cl 

m 


= -/\ 


dd>  I 


d^  ' . 

1 


Thus,  for  (j)  < tt/2,  condition  i)  implies 


= ^ f I (6')  d6'=  it  , <p  < -n/2 
2 J pv  pv  ra 


(C-3) 


where  i is  the  average  value  of  it  ..  For  <)>  > ti/2  , and  noting 

pv  ® pvi  m 

that  |d(i)/dll|,  = li  ■ / j . conditions  li)  and  iii)  imply 
' '1  ml  pvi 


!>-  = J 

c 2 


■n/4>-l 


f H (S')  d6'  + 7t  f It  (6') 

J pv  J 


-1 


Tr/<p  -1 
m 


d6’ 

pv^"  ^ <1.  (6') 

m 


or 


H = I 

c 2 


it  TT 

pv 


m 


m 


f d6’ 

J 1 + A' 


tt/(P  -1 
m 


TT  J, 


^ = O A 

c 2 1^1 


£V 


ti/2 


(C-4) 


(05) 


Since  N = S/As  = S/(t  , Equations  (C-3)  and  (C-5)  determine  N in  terms 
c c 

of  4 and  I = di  /<  I ddi/dt  1 ^ . A sketch  of  the  behavior  of  I with  f 
m pv  m ‘ ' cm 

Is  shown  in  Figure  C-2.  Only  when  is  small  is  a decorrelation  length 

based  on  the  autocorrelation  of  f the  same  as  that  based  on  the  auto- 
correlation of  Aa^^. 

Consider  the  case  where  the  variation  of  d>  due  to  surface 
roughness.  For  this  case,  d>  = 2k  h = 4tt  h /3.  Equation  (C-3) 
becomes 


I ==  I = Ip 
c pv  m 


or 
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(C-6) 


I = 
c 


h < X/8 
m 


Equation  (C-5)  becomes 


I 

c 


8 < 


^1>L 


1 + in 


h > \/8 
m 


(C-7) 


2 

and  N is  given  by  S/i  . 

c 

The  quantity  h is  related  to  the  rms  height  variation  (from 
zero)  by 


h = h /v’3  . (C-8) 

rms  m 


4 


;9. 
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Appendix  D.  THE  BEHAVIOR  OF  ^ 

Consider  (p  > tr , but  N sufficiently  large  that  .})  /N  «1.  Using 
m 

Equations  (47)  and  (48)  with  Equation  (49)  and  dividing  numerator 

2 

and  denominator  by  sin  i))^  gives 


rms  m 


(1  - COS.},^)^  (1  - cost^)^ 

2 - + 1 + ' r 47  " ' 3N 

sin  <()  16  sin  4> 

hn L m J 

r (1  - 

— y-  - 3N 
4 sin 


From  trigon-  met  r ic  identities,  1 - cos((i^  = 2 sin  ^tid 

slnit  = 2 sint  /2  cos4>  /2.  Therefore, 
m mm 


2 2 

(1  - cos()i  ) sin  <|i  /2 
m m 


sin^^ 


cos  <))  /2 
m 


which  has  a minimum  value  of  zero  at  <t>^  - 2nTr,  n 1,  2,  3,  ...) 

Hence,  using  Equations  (C— 7)  and  (C— 8)  and  the  fact  that 

p = 4 TT  h /A  , one  obtains 
m m 


(61)  2 r S/T  h 

'■'^^^rms  ^ V ^ jn_  rms  L + s,n 


<I>  3N  4S 


dH 


(D-1) 


where  X = 2/3*  h /n. 

rms 

2 2 

When  pJ2  rnr  + tt/2,  n = 0,  1,  2,  3,  ...,  then  sin  {pj2)lcos  iPjD 
-t-  00,  and 


(D-2) 


Hence,  for  <|>  > n but  ^ /N  « 1 , oscillates  between  the 

^ tn  m rni.s 

limits  given  by  Equations  (D-1)  and  (D-2). 
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Now.  for  (f  so  large  (A  so  small)  that  41  /N  » 1 , Equation  (49) 
m ^ 

hei'omes  [using  Equations  (47)  and  (48)] 

1/2 

(61)  fl  +_^sin^(4i  /2)/cos^(4>  /2)j 

^ rms  ^ L W ^ (D_ 

T sin^(4'  /2)/cos^(4>  /2) 

4 m m 

which  oscillates  rapidly  with  A and  is  bounded  by  1 and  /l7/5.  If  a 

spread  in  6 exists,  then  cos(4i  /2)  ->■  0 for  large  4 (or  small  A) 
m m m 

and  (61)  /<I^  -*■  1 with  no  oscillations, 

rms 

Whenever  A . then  4>  0 and,  hence,  ('5D  > *■  0-  Thus, 

m rms 

for  A > A (61)  /<!>  falls  off  from  unity  and  approaches  zero. 

1 rms 

The  sketch  illustrating  these  mathematical  points  in  shown  in 
Figure  4, 
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LIST  OF  SYMBOLS 


Far-field  electrl-  fit*  Id  amp  1 i t luli.* 

Variation  of  A from  tin*  averaei*  A 

a 

Electric  field  amplitude  of  liglit  incident  on  (reflected 
from)  the  constructed  surface  of  the  target 

Electric  field  amplitude  of  reflected  light  at  constructed 
surface  without  phase  factor  due  to  curvature 

Variation  of  a,  from  the  average  a, 

1 la 

Constant  for  given  wavelength  and  surface  shape  and  size 

Factor  which  depends  only  on  constructed  surface  size  and 
shape  at  any  point 

Constant-of-proportional ity  between  1 and  A*A 

Distance  between  actual  target  surface  and  the  constructed 
surface  along  the  path  of  reflected  i igiit 

Deviation  of  surface  "height"  from  the  constructed  surface 
(for  which  <ij)>  = 0) 

Maximum  value  of  the  cel 1-averaged-h  over  the  exposed 
target  surface 

Far-field  light  intensity  upon  reflection 
Variation  of  I from  the  average  <I> 

Imaginary  part  of  — 

Magnitude  of  light  wave  vector 

Average  decorrelation  length  • t. . ■ based  on  autocorrelation 
of  Aaj^  * 

"Peak-to-va  1 ley"  distance  for  Aa,  at  small  ')) 

1 m 

Number  of  decorrelation  area  cells  on  exposed  target  surface 
Surface  normal  vector 


Reflectance  of  target  surface 


iip  Path  length  difference  between  incident  and  reflected 

light  at  constructed  surface 

R Distance  from  observer  to  nearest  point  on  target 

r Distance  from  observer  to  a point  on  the  constructed 

surface 

Re  Real  part  of  — 

S Exposed  surface  area  of  target 

5 Cross-section  area  of  target 
a 

a Surface  orientation  angle  for  target 

r Ratio  (61)  /<1> 

rms 

6 Fluctuation  of  — 

6!  Range  factor  for  variation  of 

6 Angle  between  _r  and  n^  (for  large  distances) 

4i  Phase  change  upon  reflection  due  to  detailed  surface 

property  fluctuations 

0 Maximum  value  of  * over  the  exposed  target  surface 

m 

^ Angle  between  light  reflected  through  the  constructed 

surface  and  the  line  from  that  point  to  observer 
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